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ABSTRACT 

The cooling rate of young neutron stars gives direct insight into their internal makeup. Although the 
temperatures of several young neutron stars have been measured, until now a young neutron star has 
never been observed to decrease in temperature over time. We fit 9 years of archival Chandra ACIS 
spectra of the likely neutron star in the ~330 years old Cassiopeia A supernova remnant with our 
non-magnetic carbon atmosphere model. Our fits show a relative decline in the surface temperature 
by 4% (5.4cr, from 2.12 ± 0.01 x 10 6 K in 2000 to 2.04 ± 0.01 x 10 6 K in 2009) and observed flux (by 
21%). Using a simple model for neutron star cooling, we show that this temperature decline could 
indicate that the neutron star became isothermal sometime between 1965 and 1980, and constrains 
some combinations of neutrino emission mechanisms and envelope compositions. However, the neutron 
star is likely to have become isothermal soon after formation, in which case the temperature history 
suggests episodes of additional heating or more rapid cooling. Observations over the next few years 
will allow us to test possible explanations for the temperature evolution. 

Subject headings: dense matter — neutrinos — stars: neutron — stars: pulsars — supernovae: indi- 
vidual (Cassiopeia A) — X-rays: stars 



1. INTRODUCTION 

The internal composition and s tructure of neutron 
stars (NSs) remains unclear (e.g. lLattimer fe Prakashl 
I2004D . Areas of uncertainty include whether exotic con- 
densates occur in the NS core, the symmetry energy and 
thus proton fraction in the core, the behavior of super- 
fluidity among neutrons and protons, the conductivity of 
the NS crust, and the chemical composition of the outer 
envelope. NSs are heated to billions of degrees during 
supernovae, and cool via a combination of neutrino and 
photon emission. Observing the cooling rates of young 
NSs is a critical method to constrain the uncertainties 
(seelTsurutalll998trYakovlev fc Pethickll2004t IPage et al.l 
120061 for reviews). 

To date, observations of young cooling NSs have 
been restricted to measuring the temperature of indi- 
vidual NSs at one point in time. As NSs may dif- 
fer in their mass, envelope composition, etc., a mea- 
surement of the cooling rate of a young NS is needed 
to determine its cooling trajectory. Since neutrino ra- 
diation (rather than the observed photon radiation) 
is the dominant source of cooling during the first ~ 
10 5 years, measurements of cooling rates during this 
time require measuring a temperature decline over time. 
No young NS has previously been observed to cool 
steadily over time. Though the ~ 10 6 -years-old NS 
RX J0720.4— 3125 has shown temperature variations of 
~10% over sa 7 years tide Vries et al.l 120041: IHohle et al.l 
120091 ). this variation is ascribed to either a glitch-like 
event or precession of surface hot spots (Haberl et al.l 
l2lMlvaTKerkwiik et al1l2007t IHohle et al.H2009D . Mag- 
netars, such as 4U 0142+61, have shown temperature 
variations, along with changes in their pulsed fraction 



due to magnetic field reconfiguration events. 

The compact central object at the center of the Cas- 
siopeia A (Cas A) supernova remn ant was discovere d 
in Chandra's first-light observations (jTananbauml ll999') . 
and quickly identified as a likely NS, which we assume 
here. It is presently the youngest-known NS, as the 
remnant's estimated age is « 330 years (jFesen et al.1 
12000). It is r elatively close-by (d = 3.4±°;^ kpc, 
iReed et aLlll995D and the supernova remnant has been 
well-studied, with over a megasecon d of Chandra ACIS 



observations sp read o yer 10 years (IrTwang et al. 120041 : 
iDeLanev et all l200l iPatnaude fc Fesenl 120071 12009( 1. 



variations, along witn c hanges in their pulsed traction 
and pulse profile (|Dib et al.l I2007T ). but these are likely 
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However, its spectrum (modeled as a blackbody or a 
magnetic or non-magnetic hydrogen atmosphere) was 
inconsiste nt with emission f rom the full surf a ce of 
the NS (IPavlov et all f2000l iChakrabartv et al.l I200H 
iPavlov fc Lunal [2009( 1 . Timing investigations using the 
Chandra HRC and XMM-Newton have failed to iden- 
tify pulsations down t o a pulsed fraction leve l of <12% 
([Murray et al.l 120021: [Mereghetti et al.l 120021: iRansoml 
I2002t lHalpern fc GotthelJ 120101 ) . Indicating that the 
emission is probably from the entire surface. These ap- 
parently contradictory observations arc reconciled by the 
discovery that an unmagnctized (B < 10 11 G) carbon at- 
mosphere provides a good fit to the Chandra ACIS data, 
with the e mission arising from the entire surface of the 
C as A NS (IHo fc Heinkef l2009h. 

IPavlov et al.1 (12004( 1 examined two long ACIS obser- 
vations (50 ks each) of the Cas A NS from 2000 and 
2002, along with several short (2.5 ks) calibration ob- 
servations, finding no significant changes in flux. Upon 
re-examination of archival Einstein and Rosat data, 
the NS was only barely detected, and thus could no t 
be used to sea r ch fo r variability (jPavlov et al.1 12000( 1 . 
IPavlov fc Lunal (12009( 1 mention that the flux measured 
in their 2006 observation is slightly lower than reported 
previously, but do not at tempt to dete r mine whether the 
difference is real. Before IHo fc Heinkei (12009(1 . it was not 
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expected that the emission arises from the entire surface 
of the NS, so further serious searches for temperature 
variations were not undertaken. Here we utilize the full 
Chandra ACIS archive of Cas A NS observations to mea- 
sure the temperature changes from 2000 to 2009. 

2. X-RAY ANALYSIS 

We analyzed all Chandra ACIS-S exposures without 
gratings, longer than 5 ks, of Cas A, listed in Table 1. 
We also analyzed the zeroth-order grating spectrum from 
ObsID 1046, which was taken in 2001. Although the 
fit to a carbon atmosphere model spectrum was good, 
the derived temperature (logT s = 6.282^'oo4) is signifi- 
cantly lower than all other ACIS measurements. This is 
likely due to calibration differences between the zeroth- 
order grating observation and observations without grat- 
ings. Similar cross-calibration uncertainties prohibit di- 
rect comparison of Chandra HRC or XMM-Newton ob- 
servations with Chandra ACIS observations. The HRC-I 
observations lack spectral information; the HRC-I team 
also uses the Cas A NS as a quantum efficiency calibra- 
tion sourcqf], which we suspect may be negatively im- 
pacting the HRC-I calibration. XMM-Newton observa- 
tions suffer substantially increased background from the 
supernova remnant. 

ObsID 6690 was taken using a subarray mode to al- 
leviate the effects of pileup on the Cas A NS; all other 
data were taken in full-frame mode, with frame times of 
3.04 or 3.24 s. We used CIAO 4.2 (with CALDB 4.2.1) 
to reprocess the observations with current calibrations, 
extract spectra and create responses. We used a 4-pixel 
(2.37") radius region for source extraction, and an annu- 
lus from 5 to 8 pixels for bac kground. Our source re gion 
is slightly larger tha n used bv lPavlov fc Lunal (|2009f ) and 
IHo fc Heinkel (J2009D . giving a more complete flux esti- 
mate when the point-spread function is asymmetric (the 
NS was slightly off-axis in most observations) . Most data 
were taken in GRADED mode, so (apart from ObsID 
6690) we could not correct the data for charge-transfer in- 
efficiency The time-dependent ACIS quantum efficiency 
degradation is modeled, but has a small effect on our 
analysis due to the few counts below 1 kcV. 

The nature of our analysis requires that we consider 
possible instrumental effects on the effective area care- 
fully. We have identified the following possible effects: 
contaminant effects on low-energy quantum efficiency 
(QE); increasing charge-transfer inefficiency effects on 
QE; bad pixels/columns; and pileup. 

Charge-transfer inefficiency can alter the "grade" des- 
ignation of events from "good" (likely X-rays) to "bad" 
(usually cosmic rays) grades, leading to deletion of good 
events. However, the Chandra X-ray Center maintains 
accurate calibration files for the S3 chip without CTI cor- 
rection (acisD2000-01-29qeuN0005.fits), which addresses 
GRADED mode data, and in any case the QE is signifi- 
cantly affected only below 1 kc\0- A molecular contam- 
inant has been building on the ACIS detector, reducing 
the QE at 1 keV by ~10% over the Chandra mission^ 
However, the effects of this contaminant are negligible 
above 2 keV, and are now calibrated across the detector 

3 cxc.harvard.edu/ccr/proceedings/07_proc/prcscntations/possonb 

4 CXC HelpDesk ticket 12871. 

5 http://cxc.harvard.edu/cal/Acis/Cal_prods/qcDcg 



and through time; uncertainties in changes of effective 
area with time are believed to be <3% for energies above 
0.7 keV. If this affected our data, we should expect to see 
greater variations in the lower-energy than higher-energy 
data, which is not the case (see below). 

Bad pixels may affect the inferred QE by removing 
good datsfl A bad pixel region at CHIPX=495-499, and 
the node boundary at CHIPX=512-513, were crossed 
by the dither pattern of the Cas A NS in the Hwang 
datasets in April/May 2004. The CHIPX=496-498 bad 
pixels were not telemetered to the ground (bias values of 
4095), so it is not possible to check the effective area cal- 
ibration by extracting data including the bad pixels. We 
analyzed this data using the appropriate bad pixel lists, 
and though the average temperatures were within the 
range of our other data, we saw significant (~ 2%) varia- 
tions in the fitted temperatures correlated with changes 
in the sky position of the bad pixels on short (~week) 
timescales. We suspect these changes are due to the re- 
sponses incompletely adjusting for the effects of the bad 
pixels, and therefore do not include data with bad pixel 
regions crossing the Cas A NS in our study of the tem- 
perature variations (though we list the results for com- 
pleteness) . 

Pileup is the recording of two photons during one fram- 
etime as one event, leading to changes in the spectrum 
and th e rejection of some events due to a change in their 
grade (jDavisI 120011 ) . Although the absolute effects of 
pileup remain uncertain, most of our observations suf- 
fer the same level of pileup. ObsID 6690 suffers much 
less pileup due to its short (0.3 s) frame time; the first 3 
observations have a slightly longer frame time (3.24 s vs. 
3.04 s). Thus we suspect that ObsID 6690 may have a 
systematic shift in T compared to the other observations 
(though it is probably the most accurate in an absolute 
sense); we show results from ObsID 6690 but exclude it 
from fits to the temperature trends. 

After some experimenting, we chose to group 50-ks 
datasets by 200 counts, and to increase the grouping for 
longer exposures to produce a similar number of bins. 
Wc merge ObsID 9117 with 9773, and ObsID 10935 
with 12020, as the datasets arc short (^25 ks) and ad- 
jacent in tim e. We use a similar spectral fit to that in 
I Ho fc Heinkel (|2009l ). a model cont aining photoelectri c 
absorption (with abundances from iWilms et al.l 120001) , 
an unmagnetized car bon NS atmosphere , and scattering 
by interstellar dust (Prcdehl et al. 2003), all convolved 
with the Davis pileup model (jDavisI 12001 in XSPEC. 
We fix the distance to 3.4 kpc. 

There is degeneracy among the fit parameters for the 
NS mass M and radius R, interstellar absorption JVh, 
and surface temperature T s . Since M, R and Nh are not 
expected to vary between observations, these are held 
constant at the best-fit values in order to explore varia- 
tions in T s . If we exclude the bad-pixel- affected data, 
the best-fit gives M = 2.01 M , R = 8.3 km, and 
A H = 1.82 x 10 21 , while the best-fit including this data 
gives M = 1.65 M , R = 10.3 km, and N H = 1.74 x 10 21 . 
We use the best-fit values for M, R and Ah derived from 
including all the data, as this provides substantially more 
information on the spectral curvature (and a more be- 
•owfeyablc NS mass), but wc do not consider the temper- 

6 http://cxc.harvard.edu/cal/Acis/Cal_prods/badpix/index.html 
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Energy (keV) 

Fig. 1. — Illustrative spectral fit of Chandra ACIS data to our 
non-magnetic carbon atmosphere model spectrum, with all tem- 
peratures forced to be equal. Upper panel shows the data and 
model. Lower panel shows the ratio of data/model, with the dif- 
ferent datasets marked: 2000: (black) plain crosses, 2002: (red) 
asterisks, Feb. 2004: (green) filled squares, 2006: (blue) circles, 
merged 2007: (cyan) stars, merged 2009: (magenta) open squares. 
The 2006 data has a higher countrate (upper panel) due to its 
lower pileup fraction. The change of countrate from early (2000- 
2002) to later (2006-2009) spectra, and its spectral uniformity, can 
be clearly seen in the ratio plot. 



atures derived from the questionable data in our analy- 
sis of temperature variations. Using the best-fit values 
without the bad-pixcl-afFectcd data gives similar results 
to those described below, except that the temperatures 
arc all shifted slightly higher by about the same amount. 
The allowed range of absolute temperatures, and the re- 
lationships between the fit parameters are explored in 
detail in Yakovlev et al. (in prep). We al low the grade 
migration parameter a in the IDavisI (|2001[) pileup model 
to vary between observations with different frame times, 
giving values of a = 0.27 ± 0.06 for the 3.24 s frame time 
observations, a = 0.24 ± 0.05 for the 3.04 s observations, 
and a < 0.62 for ObsID 6690 (all 90% confidence). Our 
results are not substantially changed by requiring a to be 
fixed across all observations. We quote la confidence er- 
rors for T s and Lboi in Table 1, for convenience in fitting 
the temperature variations. The observed flux change 
(~21% in absorbed flux) docs not depend on the choice 
of spectral model. 

We find that requiring the same NS temperature 
for the (non-bad-pixel-affected) observations produces a 
poor x 2 (= 184 for 131 degrees of freedom) with a null hy- 
pothesis probability ofl.6xl0 -3 . Allowing the NS tem- 
perature to vary reduces the \ 2 to 113 for 126 degrees of 
freedom, which gives a probability of 79% for an accept- 
able fit. An F-test finds a statistic of 15.8 and probability 
of 4 x 10~ 12 , indicating that the additional parameters 
substantially improve the fit. Figure 1 shows the best-fit 
spectral fit with a constant temperature (excluding the 
bad-pixel- affected data). Substantial differences in the 
data/model ratio are clearly evident; these differences 
increase slightly at higher energies (as expected for tem- 
perature variations) from early to later observations. If 
the molecular contaminant or charge-transfer inefficiency 
were the primary cause of the variation, we would expect 
the data/model ratio to vary principally at low energies, 
which is not seen. 
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Fig. 2. — Surface temperature T B of the Cas A NS, obtained from 
spectral fits of Chandra observations, as a function of time (crosses 
and triangles indicate best-fit values with la errorbars) . Curves are 
fits (to well-calibrated data, the 5 red crosses; see text) with a sim- 
ple model for NS cooling by modified Urea or pair breaking and 
formation (MU/PBF; long-dashed) or condensate Urea with an 
iron envelope (CU; short-dashed) or a fully accreted light-element 
envelope (CUa; dotted) after a long delay in thermal relaxation, 
while the solid curves are for cooling by modified Urea with re- 
laxation shortly [« 100 y (upper) and m 20 y (lower)] after NS 
formation (see text). 

Figure [2] shows the temperature variation over the al- 
most ten years the Cas A NS has been observed by 
Chandra. We show the five well-calibrated tempera- 
ture measurements (as well as, for reference, the 2006 
subarray measurement and the 2004 measurements af- 
fected by bad pixels). A clear decrease of 3.6 ± 0.6% in 
T s (15 ± 4% decrease in bolometric luminosity) is evi- 
dent during this period, temperature drop of 1.5 ± 0.5% 
within only 21 days. This gives a temperature evolution 
timescale T s /(AT S /At) of ~ 280 y. 

3. DISCUSSION 

A detailed analysis of the Cas A NS temperature evolu- 
tion is beyond the scope of this paper. Work is underway 
on investigating the thermal history of the Cas A NS in 
the context of neutrino cooling calculations (Yakovlev et 
al., in prep.). Here we briefly describe NS cooling theory 
and provide a simpli fied model to compare with the ob- 
served evolution fseelTsurutalll998UYakovlev fc PethickJ 
[200l iPage et al.ll2006l for review) . 

The long-term thermal history of a NS is determined 
by the neutrino luminosity and heat capacity of the 
core and the composition (i.e., thermal conductivity) of 
the surface layers. At very early times, the core cools 
rapidly via neutrino emission while the temperature of 
the thermally-decoupled crust remains nearly constant. 
A cooling wave travels from the core to the surface, bring- 
ing the NS to a relaxed, isothermal state. Depending on 
the properties of the crust, the relaxation time can tak e 
- 10 - 100 y (jLattimer et alJll99l iGnedin et afll200l . 
For the next ~ 10 5 — 10 6 y, surface temperature changes 
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reflect changes in the interior temperature as neutrino 
emission continuously removes heat from the star. 

Let us assume that the Cas A NS has become 
thermally-relaxed and that the observed temperature 
decline is due solely to neutrino emission. The ther- 
mal evolution of a young NS is then governed by the 
thermal balance equation C(dT/dt) « — L„, where T 
is the interior temperature, C is the total heat ca- 
pacity, and L„ is the total neutrino luminosity. The 
ratio Ly/C thus determines the rate of temperature 
change. The heat capacity for a non-superfluid NS is 
C ~ 10 38 T 8 ergs K" 1 , where T 8 = T/10 8 K. For sim- 
plicity, we consider L v to be given by a single neutrino 
process, i.e., either the slow modified Urea (MU) pro- 
cess with L^ u ~ 10 32 T 8 ergs s~ x or a fast condensate 
Ure a (CU) process with . L™ ~ 10 36 T S 6 ergs s " 1 (see, 
e.g., lYakovlev fc Pethickl l200l iPage et alJl200l for re- 
view); note that the nucleon direct Urea process would 
result in a temperature below that observed within one 
year after becoming isothermal. The evolution equation 
then results in Tg(t) ~ To(nt)~ n , where t is in years, 77 
encapsulates the uncertainties in the coefficient of L v /C, 
and n MV = 1/6 and T MU = 9 and n cu = 1/4 and 
Tq V = 3. Note that we assume the current temperature 
is much lower than the initial temperature. To convert 
the evolution of the interior temperature to one for the 



surface temperature, we use the relation T s g « 1 . 1 T t 



11/20 

for an iron envelope (|Gudmundsson et al.l Il982f ) and 

T s6 p a 1.8T g for a (fully accret ed) light-element enve- 
lope (jPotekhin et al .11 19971 see also lPotekhin et al.ll2003h . 
where T s6 = T s /10 6 K. We thus obtain 



T a6 (t) = T a0 (r]ty 



(1) 



11/120 and 7$ 



MU 



where a MU 

and T s o = 2 for an iron envelope and a 

and T# Ul 



4 and a cu = 11/80 

MUa = 17/168 

17/112 and T^ Ua = 4 for a 



7 and a CUa 
fully accreted envelope. 

In Fig. [2j we show the predictions of the temperature 
decline given by cq. (JlJ, where t = (t — r x ), t is the 
year, and r x is the approximate year when the star be- 
comes isothermal (after which our simple thermal evo- 
lution scaling is valid). The deviation from unity of 
77 (which encompasses the uncertainties in our knowl- 
edge of L v /C) is a rough measure of the likelihood that 
a given process and composition is responsible for the 
temperature evolution seen in the Cas A NS. If we as- 
sume the NS took a long time to thermally relax (so 
that it only became isothermal recently), a fit to the 
observations (excluding the bad-pixcl-afiectcd data and 
the 2006 subarray mode data) for each neutrino emis- 



sion process and envelope composition yields \ 2 ~ 1-4 
for 3 degrees of freedom. For slow cooling and an ac- 
creted envelope, we find 7?(MUa)= 6000, which requires 
a significantly higher neutrino luminosity and/or lower 
heat capacity than traditionally considered; in other 
words, the observed T s requires a T that is too low 
for slow processes to achieve in 300 y. On the other 
hand, r/(MU)= 50, r/(CU)= 0.02, and 7/(CUa) = 2, and 
thus are possibile scenarios. We have not taken into 
account the large (suppression) effects of superfluidity 
on the neutrino luminosity and nucleon heat capacity 
because of the relatively-unknown critical t emperature 
at w hich these effects begin to occur (see IPage et all 
12004 and references therein). However, we note that 
if strong supcrfluid pairing exists in the NS, then neu- 
trino emission by pair breaking and formation (PBF) 
produces Ll BF > l0Lf v (|Gusakov et alJl200l : this re- 
sults in ??(PBF)< 5. We also find r x (MU/PBF)=1980, 
r x (CU)=1968, and r x (CUa)=1964, i.e., - 300 y after the 
supernova. Despite the good fit to the data, the very long 
relaxation time makes this scenario questionable. 

We also show two (MU) cooling curves in Fig. [2j 
which assume that the NS became thermally-relaxed 
early (« 20 y and « 100 y) after formation; no fit is 
done, and n w 5 is assumed. In this case, the clear 
deviation from a single neutrino cooling curve suggests 
that a transient heating (or cooling) episode occurred. 
Transient heating can originate from external causes, 
e.g., accre t ion of circumstellar gas or asteroids (e.g. 
[Jural [200l iCordes fc Shannon! 120081) . or from internal 
sources, e.g., by the same mechanism that is responsible 
for pulsar glitches (Ivan Riper et all 119911 : iSeward et al.l 
l2000HHelfand et all 1200 ll and references therein). Fur- 
ther monitoring of the Cas A NS temperature evolution 
should allow us to test the possibility of transient heating 
via a return to quiescent neutrino cooling. 
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TABLE 1 
Spectral Fits to Cas A Neutron Star 



Start time 
CUT) 


MJD 


Exposure 

ks 


Frame time 

s 


logT s 
K 


L hoh 1 eV-10 keV 
(xlO 33 ergs s _1 ) 


ObsID 


2000-01-30 10:40 


51573.74 


50.56 


Directly comparable 
3.24 


data 

fi q 9C :o+0.0021 

o.a^oo_ 0021 


7 qcc+0.15 
'• M0 -0.15 


Hi 


2002-02-06 06:22 


52311.56 


50.3 


3.24 


fi 0907+O.OO21 
o.ozo(_ 0021 


7 7 n+o.is 

'•' M -0.15 


1952 


2004-02-08 17:41 


53044.03 


50.16 


3.24 


fi q 17 n+ 00020 

o.o±(u_ 00 2o 


7 OO+0.13 
' •°°-0.13 


5196 


2007-12-05 22:00 
2007-12-08 12:34 


54441.36 

a 


25.18 
25.17 


3.04 
3.04 


fi qi 97+0-0018 

O.Ol^(_ 001g 

a 


7-05±8;i? 

a 


9117 
9773 


2009-11-02 22:24 


55138.59 


23.58 


3.04 


6 3104+ 00019 

D.J1U4_ 0Qlg 


a nn+0.12 

6.90I 011 


10935 


2009-11-03 22:32 


b 


22.68 


3.04 


6 


b 


12020 


Other data 


2004-04-14 19:47 


53110.79 


166.72 


3.04 


6 S196+ lm,ii 

D.J1JD_ 0Q11 


7 50+ U ' UY 
'• au -0.05 


4638 


2004-04-18 21:18 


53114.14 


42.84 


3.04 


fi q 1fi 7+0.0020 
O-JJ-O' -0.0021 


7 30+ ' 12 
'••ju_ 12 


5319 


2004-04-20 08:41 


53116.20 


145.38 


3.04 


6 ^164+ - 0012 
D.J1D4_ 0Q12 


7 90+O.O7 
l-^S-0.07 


4636 


2004-04-22 18:40 


53118.74 


165.66 


3.04 


6 3166+ 00011 

D.J1DD_ 0Q11 


7 oq+0.07 
'••JU-0.07 


4637 


2004-04-25 09:37 


53120.87 


80.13 


3.04 


6 S097+ - 0015 
D - JUJ ' -0.0015 


6 S^+O- 09 


4639 


2004-04-28 05:43 


53124.11 


150.59 


3.04 


6 3091 + - 0011 
D.juyi_ 0011 


R 09+0. 06 


4634 


2004-05-01 00:44 


53126.82 


136.82 


3.04 


fi S14Q+0-0° 12 
D.ji4y_ 0012 


7 1O+0.07 
'• ly -0.07 


4635 


2004-05-05 22:59 


53131.28 


55.11 


3.04 


fi o 19c; +0.0018 


7 03+ ' 11 
'• u,3 -0.11 


5320 


2006-10-19 08:18 


54027.75 


70.18 


0.34 


fi o n7q +0.0009 
D - JU ' y -0.0014 


fi 74+O.O6 
D -'^-0.08 


6690 



Note. 



N H = 1.74 x 10 21 



Carbon atmosphere spectral fits, with fixed distance of 3.4 kpc, neutron star mass of 1.648 A/q, radius of 10.3 km, and 



(best-fit). Temperature and luminosity errors arc lcr confidence for a single parameter. MJD values are for the 



midpoints of the observations, or the weighted midpoints of merged sets. 

a Combined with ObsID 9117; values for MJD, logT, and Lbol are f° r the merged set. 

logT, and I/bol are f° r the merged set. 



Combined with ObsID 10935; values for MJD, 
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